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I. INTRODUCTION 



The importance of integrable nonlinear partial differential equations (PDEs) in 1+1 dimensions 
in applications to nonlinear physics can hardly be overestimated. Their importance partially stems 
from the fact that it is always possible to obtain certain explicit solutions, called solitons, by some 
algebraic procedure. At present, there is a wide range of literature concerning integrable nonlinear 
PDEs and their soliton solutions (see, for instance, Refs. 0,0,0,0] and the references therein). The 
reader familiar with the inverse scattering transform method knows that it is zeros of the Riemann- 
Hilbert problem (or poles of the reflection coefficients in the previous nomenclature) that give rise 
to the soliton solutions. These solutions are usually derived by using one of the several well-known 
techniques, such as the dressing method [ll 0], the Riemann-Hilbert problem approach 0, 0], 
and the Hirota method (see [l]). In the first two methods, the pure soliton solution is obtained by 
considering the asymptotic form of a rational matrix function of the spectral parameter, called the 
soliton matrix in the following. It is known,that t he g eneric case of zeros of the matrix Riemann- 
Hilbert problem is the case of simple zeros 0, 0, lid Uli 13 ] (see also Ref. 13|). A single simple 
zero produces a one-soliton solution. Several distinct zeros will produce multi-soliton solutions, which 
describe the interaction (scattering) of individual solitons. As far as the generic case is concerned, 
there is no problem in the derivation of the corresponding soliton solutions. 

However, in the non-generic cases, when at least one higher-order (i.e. multiple) zero is present in 
the Riemann-Hilbert problem, the situation is not so definite. Higher-order zeros must be considered 
separately, as, in general, the soliton solutions which correspond to such zeros cannot be derived from 
the known generic multi-soliton solutions by coalescing some of the distinct simple zeros. This is clear 
from the fact that a higher-order zero generally corresponds to a higher-order pole in the soliton matrix 
(or its inverse), which cannot be obtained in a regular way by coalescing simple poles in the generic 
multi-soliton matrix. The procedure of coalescing several distinct simple zeros produces only higher- 
order zeros with equal algebraic and geometric multiplicities (the geometric multiplicity is defined as 
the dimension of the kernel of the soliton matrix evaluated at the zero), which is just the trivial case 
of higher-order zeros. For instance, if the algebraic multiplicity is equal or greater than the matrix 
dimension, then such coalescing will produce a higher-order zero with the geometric multiplicity no 
less than the matrix dimension, which could only correspond to the zero solution instead of solitons. 
Thus the soliton matrices corresponding to the higher-order zeros of the Riemann-Hilbert problem 
require a separate consideration. 

Soliton solutions corresponding to higher-order zeros have been investigated in the literature before, 
mainly for the 2 x 2-dimensional spectral problem. A soliton solution to the nonlinear Schrodinger 
(NLS) equation corresponding to a double zero was first given in Ref. 1J] but without much analysis. 
The double- and triple-zero soliton solutions to the KdV equation were examined in Ref. 0] and 
the gene ral multiple-zero soliton solution to the sine-Gordon equation was extensively studied in 
Ref. 116] using the associated Gelfand-Levitan-Marchenko equation. In Refs. , E3 , higher-order 
soliton solutions to the NLS equation were studied by employing the dressing method. In |19|, 
higher order solitons in the Kadomtsev-Petviashvili I equation were derived by the direct method 
and the inverse scattering method. Finally, in our previous publication 22] we have derived soliton 
matrices corresponding to a single elementary higher-order zero — a zero which has the geometric 
multiplicity equal to 1. Our studies give the general higher-order soliton solutions for the integrable 
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PDEs associated with the 2x2 matrix Riemann-Hilbert problem with a single higher-order zero. 
Indeed, any zero of the 2 x 2-dimensional Riemann-Hilbert problem is elementary since a nonzero 
2x2 matrix can have only one vector in its kernel. 

However, the previous investigations left some of the key questions unanswered. For instance, the 
general soliton matrix corresponding to a single non-elementary zero remained unknown. Such zeros 
arise when the matrix dimension of the Riemann-Hilbert problem is greater than 2. Naturally then, 
the ultimate question — the most general soliton matrices corresponding to an arbitrary number of 
higher-order zeros in the general NxN Riemann-Hilbert problem, was not addressed. Because of these 



utions to PDEs integrable through the 
23|, the three- wave interaction system 



unresolved issues, the most general soliton and multi-soliton so 
NxN Riemann-Hilbert problem (such as the NLS equation 
I2H24, 25, 26, 27 1, and the Manakov equations 28 1) have not been derived yet. 

In this paper we derive the complete solution to the problem of soliton matrices corresponding 
to an arbitrary number of higher-order normal zeros for the general NxN matrix Riemann-Hilbert 
problem. These normal zeros are defined in Definition^ and are non-elementary in general. They 
include almost all physically important integrable PDEs where the involution property [see Eq. (@J)] 
holds. The corresponding soliton solutions can be termed as the higher-order multi-solitons, to reflect 
the fact that these solutions do not belong to the class of the previous generic multi-soliton solutions. 
Our results give a complete classification of all possible soliton solutions in the integrable PDEs 
associated with the NxN Riemann-Hilbert problem. In other words, our soliton matrices contain 
the most general forms of reflection-less (soliton) potentials in the iV-dimensional Zakharov-Shabat 
spectral operator. For these general soliton potentials, the corresponding discrete and continuous 
eigenfunctions of the iV-dimensional Zakharov-Shabat operator naturally follow from our soliton 
matrices. As an example, we consider the three-wave interaction system, and derive single-soliton 
solutions corresponding to a non-elementary zero, and higher-order two-soliton solutions. These 
solutions generate many new processes such as the simultaneous breakup of a pumping wave (^3) 
into the other two components (ui and U2) and merger of Ui and U2 waves into the pum ping w 3 wave, 
i.e., u\ + 112 + u% <-> u\ + U2 + «3- They also reproduce previous solitons in l27j ] as special 

cases. 

The paper is organized as follows. A summary on the Riemann-Hilbert problem is placed in section 
im Section 11111 is the central section of the paper. There we present the theory of soliton matrices 
corresponding to several higher-order zeros under the assumption that these zeros are normal (see 
Definition Q), which include the physically important cases with the involution property [see Eq. P|)]. 
Applications of these general results to the three-wave interaction system are contained in Section 
II VI Finally, in the appendix we briefly treat the more general case where the zeros are abnormal. 



II. THE RIEMANN-HILBERT PROBLEM APPROACH 



The integrable nonlinear PDEs in 1+1 dimensions are associated with the matrix Riemann-Hilbert 



problem (consult, for instance, Refs. Q S, S H, 0, B 0, H 0, Q H H H H E3). The 
matrix Riemann-Hilbert problem (below we work in the space of N x N matrices) is the problem of 
finding the holomorphic factorization, denoted below by $+(/c) and $I 1 (/c), in the complex plane of 
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a nondegenerate matrix function G(k) given on an oriented curve 7: 

$Z\k, x, t)$+(k, x, t) = G{k, x, t) = E{k, x, t)G{k, 0, 0)E~ 1 {k, x,t), k e 7, (1) 

where 

E(k, x, t) = exp [-A(k)x - tt(k)t). 

Here the matrix functions <&+(k) and are holomorphic in the two complementary domains of 

the complex fc-plane: C + to the left and C_ to the right from the curve 7, respectively. The matrices 
A(fc) and fl(k) are called the dispersion laws. Usually the dispersion laws commute with each other, 
e.g., given by diagonal matrices. We will consider this case (precisely in this case E(k, x, t) is given by 
the above formula). The Riemann-Hilbert problem requires an appropriate normalization condition. 
Usually the curve 7 contains the infinite point k = 00 of the complex plane and the normalization 
condition is formulated as 

$±(k,x,t) — > I, as k^oo. (2) 

This normalization condition is called the canonical normalization. Setting the normalization condi- 
tion to an arbitrary nondegenerate matrix function S(x, t) leads to the gauge equivalent integrable 
nonlinear PDE, e.g., the Landau-Lifshitz equation in the case of the NLS equation Obviously, 
the new solution $±(/c, x, t) to the Riemann-Hilbert problem, normalized to S(x,t), is related to the 
canonical solution by the following transformation 

$±(k, x, t) = S(x, t)$(k, x, t). (3) 

Thus, without any loss of generality, we confine ourselves to the Riemann-Hilbert problem under the 
canonical normalization. 

For physically applicable nonlinear PDEs the Riemann-Hilbert problem possesses the involution 
properties, which reduce the number of the dependent variables (complex fields). The following 
involution property of the Riemann-Hilbert problem is the most common in applications 

& + (k) = ^(k), k = k*. (4) 

Here the superscript "f" represents the Hermitian conjugate, and "*" the complex conjugate. Exam- 
ples include the NLS equation, the Manakov equations, and the N-wave system. The analysis in this 
article includes this involution (@J as a special case. In this case, the overline of a quantity represents 
its Hermitian conjugation in the case of vectors and matrices and the complex conjugation in the 
case of scalar quantities. In other cases, the original and overlined quantities may not be related. 

To solve the Cauchy problem for the integrable nonlinear PDE posed on the whole axis x, one 
usually constructs the associated Riemann-Hilbert problem starting with the linear spectral equation 

d x $(k, x, t) = $(k, x, t)A(k) + U(k, x, t)$(k, x, t), (5) 
whereas the t-dependence is given by a similar equation 

d t $(h, x, t) = x, t)tt(k) + V(k, x, t)$(k, x, t). (6) 

The nonlinear integrable PDE corresponds to the compatibility condition of the system (0) and (JHJ): 

d t U - d x V + [U, V] = 0. (7) 
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The essence of the approach based on the Riemann-Hilbert problem lies in the fact that the 
evolution governed by the complicated nonlinear PDE (JJJ) is mapped to the evolution of the spectral 
data given by simpler equations such as (JIJ and ()2()aj) - (J20bj) . When the spectral data is known, 
the matrices U(k, x, t) and V(k, x, t) describing the evolution of $± can then be retrieved from 
the Riemann-Hilbert problem. In our case, the potentials U(k,x,t) and V(k,x,t) are completely 
determined by the (diagonal) dispersion laws A(k) and Q(k) and the Riemann-Hilbert solution $ = 
$±(k,x,t). Indeed, let us assume that the dispersion laws are polynomial functions, i.e., 



Ji -h 

A(fc) = A o k ^ Q ( k ) = Yl B ^ j - 

3=0 3=0 



(8) 



Then using similar arguments as in Ref. |32| we get: 

u = -v{$a$~ 1 }, v = -v{$n$- 1 }. (9) 

Here the matrix function $(/c) is expanded into the asymptotic series, 

= / + £T 1 $ (1) + Ar 2 $ (2) + ..., k^oo, 

and the operator V cuts out the polynomial asymptotics of its argument as k — > oo. An important 
property of matrices U and V is that 



TrU(k, x, t) = — TrA(fc), TrV(k, x, t) = -TrQ(k), 



(10) 



which evidently follows from equation Q. This property guarantees that the Riemann-Hilbert zeros 
are (x, £) -independent. 

Let us consider as an example the physically relevant three- wave interaction system [3. 124 . 25. 27 . 
Set N = 3, 



^ Oi \ 



A(fc) = ikA, A 



a 2 
a 3 



( h o \ 



fl(Jfe) = ikB, B 



ft 2 o 
6 3 



where aj and bj are real with the elements of A being ordered: a\ > a 2 > a 3 . From equation we 
get 

U = -A(Jfe)+i[A, V = -fi(ife) + (12) 



Setting 

assuming the involution (@J), and using equation (|T2*|) in (0) we get the three- wave system: 



«i = y/ai ~ 02*12 , M 2 = a/o2 - 03*23, u 3 = V01 - 03*^3, (13) 



(9 t Mi + v\d x ui + ieu 2 u 3 = 0, 
<%m 2 + v 2 d x U2 + i£«i«3 = 0, 
9tM 3 + v 3 d x u 3 + ieu\u 2 = 0. 



(14a) 
(14b) 
(14c) 
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Here 

b 2 -bt b 3 -b 2 h-bi 
vi = , v 2 = , v 3 = , (15) 

d\ — (X2 ^2 — ^3 ^1 — ^3 

£ _ aih - a-ibi + a 2 b 3 - a 3 b 2 + a 3 bi - aifr 3 
[(oti - a 2 )(a 2 - a 3 )(ai - as)] 1 /* 

The group velocities satisfy the following condition 

= _a 1 -a 2 < Q 
Di - v 3 a 2 - a 3 

The three-wave system f!14|) can be interpreted physically. It describes the interaction of three wave 
packets with complex envelopes u\, u 2 and u 3 in a medium with quadratic nonlinearity. 

In general, the Riemann-Hilbert problem ©-(0) has multiple solutions. Different solutions are 
related to each other by the rational matrix functions T(k) (which also depend on the variables x and 

t) BSiBQ: 

$±(k, x, t) = $±(k, x, t)T(k, x, t). (18) 

The rational matrix T(k) must satisfy the canonical normalization condition: T(k) — > I for k — > oo 
and must have poles only in C_ (the inverse function T~ 1 (k) then has poles in C + only). Such a 
rational matrix T(k) will be called the soliton matrix below, since it gives the soliton part of the 
solution to the integrable nonlinear PDE. 

To specify a unique solution to the Riemann-Hilbert problem the set of the Riemann-Hilbert 
data must be given. These data are also called the spectral data. The full set of the spectral data 
comprises the matrix G(k,x,t) on the right-hand side of equation (JTJ) and the appropriate discrete 
data related to the zeros of det $+(£;) and det$I 1 (fc). In the case of involution ®, the zeros of 
det and det $z\k) appear in complex conjugate pairs, kj = k*. It is known 

010, HQQ 

(see also Ref. [l3|) that in the generic case the spectral data include simple (distinct) zeros ki, . . . ,k n 
of det$ + (/c) and ki,...,k n of det$I 1 (/c), in their holomorphicity domains, and the null vectors 
\vi), . . . , \v n ) and (vi\, . . . , (v n \ from the respective kernels: 

$+{k j )\v j ) = 0, (v j \$Z 1 (k j ) = 0. (19) 

Using the property (|10|) one can verify that the zeros do not depend on the variables x and t. The 
(x, ^-dependence of the null vectors can be easily derived by differentiation of (JT§j) and use of the 
linear spectral equations This dependence reads: 

\vj) = exp {— A(kj)x — Q(kj)t}\voj) , (20a) 
(vj\ = (v oj \ exp{A(kj)x + tt(kj)t}, (20b) 

where \voj) and (voj\ are constant vectors. 

The vectors in equations ()20a|) - (j20b|) together with the zeros constitute the full set of the generic 
discrete data necessary to specify the soliton matrix T(k,x,t) and, hence, unique solution to the 
Riemann-Hilbert problem (|T|)-(p|). Indeed, by constructing the soliton matrix T(k) such that the 
following matrix functions 

Mk) = MAOr-W, <t>'\k) = m^ik) (2i) 
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are nondegenerate and holomorphic in the domains C+ and C_ , respectively, we reduce the Riemann- 
Hilbert problem with zeros to another one without zeros and hence uniquely solvable (for details see, 
for instance, Refs. Below by matrix T(k) we will imply the matrix from equation (|21j) 

which reduces the Riemann-Hilbert problem ©-(0) to the one without zeros. The corresponding 
solution to the integrable PDE (JJJ) is obtained by using the asymptotic expansion of the matrix $(&;) 
as k — > oo in the linear equation (jSj). In the iV-wave interaction model it is given by formula (jl2j) . 
The pure soliton solutions are obtained by using the rational matrix $ = T(k). 

The above set of discrete spectral data (|19|) holds only for the generic case where zeros of det 
and det$I 1 (fc) are simple. If these zeros are higher-order rather than simple, what the discrete 
spectral data should be and how they evolve with x and t is still unknown yet. We have stressed in 
Sec. 1 that the case of higher-order zeros can not be treated by coalescing simple zeros, thus is highly 
non-trivial. In this paper, this problem will be resolved completely. 



III. SOLITON MATRICES FOR GENERAL HIGHER-ORDER ZEROS 

In this section we derive the soliton matrices for an arbitrary matrix dimension N and an arbitrary 
number of higher-order zeros under the assumption that these zeros are normal (see Definition 0). 
Normal higher-order zeros are most common in practice. In general, they are non-elementary. Our 
approach is based on a generalization of the idea in our previous paper 2^ . 



A. Product representation of soliton matrices 

Our starting point to tackle this problem is to derive a product representation for soliton matrices. 
This product representation is not convenient for obtaining soliton solutions, but it will lead to the 
summation representation of soliton matrices, which are very useful. 

In treating the soliton matrix as a product of constituent matrices (called elementary matrices in 
Ref . , see formulae (J24j) and (J27j) below) one can consider each zero of the Riemann-Hilbert problem 
separately. For instance, consider a pair of zeros k\ and k ly respectively, of $+(&) and $I 1 (/c) from 
Eq. each having order m: 

det^+(k) = (k-k 1 ) m ^(k), detfclHfc) = {k-kx) m Tp{k), (22) 

where <p(ki) 7^ and Tp{ki) ^ 0. The geometric multiplicity of k\ ik\) is defined as the number 
of independent vectors in the kernel of $ + (/ci) ($I 1 (/ci)), see (IT9|) . In other words, the geometric 
multiplicity of k\ ik\) is the dimension of the kernel space of $ + (fci) ($I 1 (/c 1 )). It can be easily shown 
that the order of a zero is always greater or equal to its geometric multiplicity. It is also obvious that 
the geometric multiplicity of a zero is less than the matrix dimension. Let us recall how the soliton 
matrices are usually constructed (see, for instance, Refs. 0,E3|)- Starting from the solution §±(k) 
to the Riemann-Hilbert problem Q-©, one looks for the independent vectors in the kernels of the 
matrices $+(^i) and $I 1 (fci). Assuming that the geometric multiplicities of k\ and k\ are the same 
and equal to r 1; then we have 

M*i)ki> = o> (vul^ih) = o, i = i,...,n. (23) 
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Next, one constructs the constituent matrix 



xi(k) = i - fr|^ ^ 

where 

J-L 

p i = ki)(^~%'(%l> K ij = (va\vji)- (25) 

Here Pi is a projector matrix, i.e., Pf = Pi. It can be shown that detxi = (k — ki) ri /(k — &i) ri [note 
that the geometric multiplicity ri is equal to rankPj. If r\ < m then one considers the new matrix 
functions 

$+(£;) = ^ + (k)xi 1 (k), $z\k) = X i(k)^z\k). 

By virtue of equations (J23|) . the matrices $+(&) and $l]j7c) are also holomorphic in the respective 
half planes of the complex plane (see Lemma 1 in Ref. 2^1 ). in addition, k\ and k\ are still zeros 



of det<3> + (/c) and det$I 1 (/c). Assuming that the geometric multiplicities of zeros k\ and k\ in new 
matrices &+(k) and $I 1 (/c) are still the same and equal to r 2 , then the above steps can be repeated, 
and we can define matrix X2{k) analogous to Eq. (ffilj) . In general, if the geometric multiplicities of 
zeros k\ and k\ in matrices 

$+(*) = ^(^xr 1 ^) • • • X7\(k), ®-J(k) = Xi-i(k) ■ ■ ■ Xi(k)^Z\k) (26) 

are the same and given by r\ (I = 1,2,...), then we can define a matrix xi similar to Eqs. (|24j) and 
(|23j) but the independent vectors \v a ) and (vu\ (i = l,...,rj) are from the kernels of $ + (fci) and 
§Z x (ki) in Eq. $Z$. When this process is finished, one would get the constituent matrices Xi(k), 
. . . , Xr(k) such that r\ +r 2 + . . . +r n = m, and the product representation of the soliton matrix T(k), 



T(k) = xn(k)--- X 2(k)xi(k), (27) 

This product representation ()27|) is our starting point of this paper. In arriving at this represen- 
tation, our assumptions are that the zeros k\ and k\ have the same algebraic multiplicity [see Eq. 
(122}) ] . and their geometric multiplicities in matrices $+(&) and $I 1 (/c) of Eq. (|26j) are also the same 
for all /'s. For convenience, we introduce the following definition. 

Definition 1 A pair of zeros k\ and k\ in the matrix Riemann-Hilbert problem are called normal 
zeros if they have the same algebraic multiplicity, and their geometric multiplicities in matrices &+(k) 
and $I 1 (A;) of Eq. Wty are also the same for all I 's. 

In the text of this paper, we only consider normal zeros of the matrix Riemann-Hilbert problem. 
The case of abnormal zeros will be briefly discussed in the Appendix. 

Remark 1 Under the involution property (@J), all zeros are normal. Thus, our results for normal 
zeros cover almost all the physically important integrable PDEs. 

Remark 2 Normal zeros include the elementary zeros of Ji3] as special cases, but they are non- 
elementary in general. 
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It is an important fact (see Ref. 22j, Lemma 2) that the sequence of ranks of the projectors p in 
the matrix T(k) given by Eq. (|2*7jl . i.e. built in the described way, is non-increasing: 



rankP n < rankP r _i < . . . < rank Pi, (28) 

i.e., r n < r n _i < . . . < r\. This result allows one to classify all possible occurrences of a higher-order 
zero of the Riemann-Hilbert problem for an arbitrary matrix dimension N. In general, for zeros of 
the same order, different sequences of ranks in Eq. (|28|) give different classes of higher-order soliton 
solutions. In Ref. 2^1 we constructed the soliton matrices for the simplest sequence of ranks, i.e., 



1,...,1. Such zeros are called "elementary". If the matrix dimension N = 2 (as for the nonlinear 
Schrodinger equation), then all higher-order zeros are elementary since rankPx is always equal to 1. 

To obtain the product representation for soliton matrices corresponding to several higher- 
order normal zeros one can multiply the matrices of the type (|2*7|) for each zero, i.e. r(/c) = 
Fi(k)T2(k) • • - Tff z (k), where Nz is the number of distinct zeros and each Tj(k) has the form given 
by formula (}2*T|) with n substituted by some rij. 

The product representation ()27|) of the soliton matrices is difficult to use for actual calculations 
of the soliton solutions. Indeed, though the representation (|2*7|) seems to be simple, derivation of the 
(x, independence of the involved vectors (except for the vectors in the first projector Pi) requires 
solving matrix equations with (x, independent coefficients. One would like to have a more convenient 
representation, where all the involved vectors have explicit (x, independence. Below we derive such 
a representation for soliton matrices corresponding to an arbitrary number of higher-order normal 
zeros. 

For the sake of clarity, we consider first the case of a single pair of higher-order zeros, followed by 
the most general case of several distinct pairs of higher-order zeros. 



B. Soliton matrices for a single pair of zeros 

Let us introduce a definition. 

Definition 2 For soliton matrices having a single pair of higher- order normal zeros (ki, ki), suppose 
r(k) is constructed judiciously as in Eq. \2T\) , with ranks r$ of matrices Pj(l < j ' < ri) satisfying 
inequality i.e., 

r n < r n -i < ■ ■ ■ < n. 
Then a new sequence of positive integers 

81 > S 2 > ■ ■ ■ > 

are defined as follows: 

s v = the index of the last positive integer in the array [ri + 1 — v, r-i + 1 — v, . . . , r n + 1 — v]. 

We call the sequence of integers {r n ,r n _i, . . . , ri} the rank sequence associated with the pair of zeros 
(ki, ki), and the new sequence {s x , s 2 , . . . , s ri } the block sequence associated with this pair of zeros. 
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Remark It is easy to see that the sum of the block sequence is equal to the sum of all ranks, 



n, 



with the latter being equal to the algebraic order of the Riemann-Hilbert zeros (ki, ki). 

For example, if the rank sequence is {3} [only one constituent matrix in (J27j) - trivial higher- 
order zero], then the block sequence is {1,1,1}; if the rank sequence is {1,1,1,1} (an elementary 
zero), then the block sequence is {4}; if the rank sequence is {2, 3, 5, 7}, then the block sequence is 
{4,4,3,2,2,1,1}. 

With these definitions the most general soliton matrices T(k) and T -1 (A;) for a single pair of higher- 
order normal zeros (k\, ki) are given as follows. This result is a generalization of our previous result 
to non-elementary higher-order zeros. 



Lemma 1 Consider a single pair of higher-order normal zeros (ki, k{) in the Riemann-Hilbert prob- 
lem. Suppose their geometric multiplicity is r\, and their block sequence is {si, 82, ■ ■ ■ ,s ri }. Then the 
soliton matrices T(k) and T^ 1 (k) can be written in the following summation forms: 



Here S v and S v are the following block matrices, 



5 * = ££ 



1 i_(i/) w _(w) 

'<? ){p; 



1+1-j 1 



{\t}),---MV))D v {k) 



l^ik-k^-i 



{\p^),... : \p^))D v {k) 



1=1 j=i ^ 11 



-1 



1 <#! ^ 



1 Hi 



D v {k) and D v {k) are the triangular Toeplitz matrices with poles 

... \ 



/ 1 

1 



D v (k) 





1 



(fc-fci) 2 (fc-fci) 





-0= L- 1 

(fc-fci) 2 (k-ki) / 



( _L 1_ 

(fc-fci) (fc-fcj) 2 



, D v (k) 



1 







1 



(fc-fci) s ^ 



1 



(fc-fci) (fc-fci) 2 



V 







1 



^-fci) / 



(29) 



(30a) 



(30b) 



(31) 



The vectors \pf^ ),{pi |, m ) {i = here are independent of k, and each of the two 



sets of vectors {\p^), Ip^ 1 ')} an d {(P^i'li •••> (P^'l} are linearly independent 
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Remark 1 If r\ — 1, the zeros ki and k\ are elementary • In this case, the above soliton matrices 
reduce to those in [22]. 

Remark 2 The total number of all |p)-vectors or (p|-vectors from all //-blocks are equal to the 
algebraic order of the zeros k\ and k\. 

Proof The representation (|29jl can be proved by induction. Consider, for instance, the formula for 
T(k). Obviously, this formula is valid for n — 1 in Eq. (|27|). where T(k) contains only a single matrix 
Xi{k)- Now, suppose that this formula is valid for n > 1. We need to show that it is valid for n + 1 as 
well. Indeed, denote the soliton matrices for n and n + 1 by T(k) and T(k) respectively, the rightmost 
multiplier in T(k) being x(k). Then we have 

T(k) = r(k)m = (l + T^ + 7rA^ + --- + 7r^) (l+ R 



k — ki (k - ki) 2 (Jfc-Jfei)V V k — k\ 

= I + + A - +...+ An - +1 , (32) 

k - ki (k-h) 2 (k - ki) n+1 V ; 

where 

f 

R=(k 1 -k 1 )P = ^2\u l )(u l \. (33) 
l=i 

Here we have normalized the vectors \ui) and {u{\ such that 

(ui\ui) = (k 1 - fci)5 M , (34) 

and f = rankR. In view of Eq. f!28|) . we know that f > ri, where r\ is the geometric multiplicity of 
k\ and fci in the soliton matrices T(k) and r _1 (A;). The coefficients at the poles in T(k) are given by 

Ai = A x + R, Aj = Aj + Aj-xR, j = 2, n, A n+1 = A n R. (35) 

Consider first the coefficients A 2 to A n+ \. The explicit form of coefficients Aj can be obtained from 
Eqs. flU, (HI, and (JS2) as 

A i = £4 ) = t'"f,'\tW^ i -,\, (36) 
I/=l I/=l 1=1 

where the inner sum is zero if s v + 1 — j < 0. Substituting this expression into (|33|) and defining the 
following new vectors in each block 

g; ] \ = (p[ u) \R, $ ) \ = (p { ; ) \r+&;}i\, j = 2, (37) 

(for blocks of size 1, s u = 1, the second formula in (pT7j) is dropped), we then put the coefficients 
A 2 , . . . , A n+ i into the required form: 

A J = X) X fo" >(Ps*+2-j-il> J = 2, . . . , n + 1, 

v=l 1=1 
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where 



kz M ) = \Qi V) ), l = l,...,s„-l, 



and s u = s u + 1, i.e., the size of each z/-block grows by one as we multiply by in formula (|32j) . 

Next, we consider the coefficient A\. Defining the vector (p^ \ = (p^J\ and utilizing the definition 
(jHZj), we can rewrite A± as 

a? = E " E itW^JR- (38) 

1=1 1=2 

To put A\ — A\ + R into the required form 

*i= £ i^'xp^i+EEi?!"')^!-,!. (39) 

/i=ri+l ^=1 Z=l 

we must define exactly one new vector \~q~^) for each z/-block [in the second term of Eq. (J39j) ] and 

r — r\ new blocks of size 1 containing 2(f — n) new vectors ) and (p^* |. Due to formulae f)35|) 
and ((SSI ) the new vectors to be defined must satisfy the following equation 

£ fi?">flf"i + £ ^?><pf"ifl = «-££ i?!"')^-^. (40) 

/i=ri+l y=l i/=l Z=2 

where the (p[ ^| definition in Eq. (|37|) has been utilized. Substituting the expression for i? into 
the above equation, we get 

E i^^^EifiXoii. (4i) 

n=n+i i=i 

where 

To show that Eq. (|4*H) is solvable, we need to use an important fact, i.e., the matrix 
M = (M Vi i), = (pJ'Vj), i/= l,...,r ls Z = l,...,fi, 

has rank ri. This fact can be proved by contradiction as follows. 

Suppose the matrix M. has rank less than r±. Then its r\ rows are linearly dependent. Thus, there 
are such scalars Ci, C2, • • • , C ri , not equal to zero simultaneously, that the vector 

n 



^=1 



is orthogonal to all |tij)'s, i.e., 

(rj\ui) = 0, 1 < Z < f. (42) 
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According to our induction assumption that soliton matrices for n have the form (|29p. we can easily 
show from the identity T(k)T~ 1 (k) = / that (p\ v ^\T~Hk-\)=0 for all 1 < v < t\ (see Thus 
(?7 1 r — 1 (fci ) = as well. According to Lemma 1 in (22J, if (r/\ is in the kernel of T~ 1 (ki) and is 
orthogonal to all |«z)'s, then (77] is in the kernel of r _1 (A; 1 ) as well, i.e., (?7|r _1 (A; 1 ) = 0. But according 
to our construction of soliton matrices [see Eq. (|2*7j) ]. the vectors (ui\ (I — 1, f) are all the linearly 
independent vectors in the kernel of T~ 1 (ki). Thus (rj\ must be a linear combination of (wz|'s. Then 
in view of Eqs. (J34"|) and (j42)l . we find that (77 1 = 0, which leads to a contradiction. 

Now that the matrix M. has rank r*i, then we are able to select vectors \q^) [y = 1, • • • , n.) such 

that ri of the f vectors are zero. With this choice of \q Sv )'s, the r.h.s. of Eq. (jH| becomes f — ri 
blocks of size 1. Assigning these blocks to the l.h.s. of (j41j) . then Eq. (}4"T]) can be solved. Hence we 
can put the coefficient A\ in the required form ([39)1 . 

Next we prove that all vectors (p-^ | (1 < v < f) in the matrix are linearly independent. 
These vectors were defined in the above proof as 

f 

= (p[^\R = j^Wfal, 1 < f < ri, (43) 
i=i 

and (p-^ | for ri + 1 < v < f are simply equal to r — r\ of the vectors u\ depending on what t\ x r\ 
submatrix of Ai has rank r±. To be definite, let us suppose the first r\ columns of the matrix A4 
have rank r\ (i.e., linearly independent). Then according to the above proof, we can uniquely select 

vectors \q s ) {y = 1, . . . , ri) such that |^) = for 1 < I < n. Thus, 

(Pil = (^l> r! + l<z/<f. (44) 

Recalling that vectors (m^I (1 < z/ < f) in the projector i? f!33|) are linearly independent, and the first 

j~i columns of matrix A^l have rank ri, we easily see that vectors (Px | (1 < ^ < r) as defined in Eqs. 
(}4"3*|) and f)44j) are linearly independent. 

Lastly, we prove that the sizes of blocks in representations ()29|) are given by the block sequence 
defined in Definition |21 An equivalent statement is that the numbers of matrix blocks with sizes 
[1, 2, 3, n] are given by the pair-wise differences in the sequence of ranks: [r x — r 2 , r 2 — r 3 , r n _i — 
r n ,r n ], where the last number in the sequence defines the number of blocks of size n. This can be 
easily proven by the induction argument using the fact that the number of new blocks of size 1 in A\ 
()35|) is given by f — n, while the sizes of old blocks grow by 1 in each multiplication as in formula 

022). 

Using similar arguments, we can prove that the representation (J29)) for T -1 (/c) is valid, and vectors 
Ip^), \pi^) are linearly independent. This concludes the proof of LemmaQ] Q.E.D. 

C. Soliton matrices for several pairs of zeros 

Next, we extend the above results to the most general case of several pairs of higher-order normal 
zeros {(ki, ki), . . . , (&jv z , k^ z )}. In this general case, the soliton matrix T{k) can be constructed as 
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a product of soliton matrices (|2"7)l for each zero by the procedure layed out in the beginning of this 
section [see Eqs. to (|2"7|)]. Thus, T(k) can be represented as 



r(k)=T 1 (k)-T 2 (k)---T Nz (k). 



(45) 



For each pair of zeros (k n ,k n ), we can define its rank sequence and block sequence by Definition^ 
either from T(k) directly or from the individual matrix T n (k) associated with this zero. It is easy to 
see that using T(k) or T n (k) gives identical results. The inverse matrix r _1 (fc) can be represented in 
a similar way. 

The product representation (J45|) for T(k) and its counterpart for T~ 1 (k) are not convenient for 
deriving soliton solutions. Their summation representations such as Eq. (|29|) are needed. It turns 
out that r(fc) and r _1 (fc) in the general case are given simply by sums of all the blocks from all pairs 
of zeros plus the unit matrix. Let us formulate this result in the next lemma. 



Lemma 2 Consider several pairs of higher-order normal zeros {(k\, ki), . . . , (k^ z , k^ z )} in the 
Riemann- Hilbert problem. Denote the geometric multiplicity of zeros (k n} k n ) as r^\ and their block 
sequence as {s^ , s% > • • • , s (1 < n < N z ). Then the soliton matrices T(k) and T~ 1 (k) can be 

written in the following summation forms: 



N Z rW N Z 

a # a 'v 

n=l v=l 



Here Sjfi and are the following block matrices, 



n=l u=l 



f(p[^\\ 



(l*W>, 



,\q { rW:\k) 



/ 



si n) = {\p ( r ] ),...Ap { ^))D^\k) 



» 



\Wi '1/ 



D^\k) and D V J' [k) are the triangular Toeplitz matrices with poles: 



D { :\k) 



i 



{k k n ) 

i 



o 
i 



\ 



{h kn)^ (h kn) 



— = — W 





1 1_ 

fcn)^ (k fert) / 



Di n \k) 



i 







V 



(46) 



(47a) 



(47b) 



WT 



( k kn ) (k kn ) ^ 
(/c k n ) 



/ 



(48) 

Vectors |p^' n ' ) ), (p£ ], |, \q\ ) (i — 1, . . . ,s^) are independent of k. In addition, for each n, 
vectors {\p[' n ^) , ...^pf 1 ' n ^)} and {(p{ ,n ^\, {pf 1 |} are linearly independent respectively. 
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Remark When there is only a single pair of zeros (ki, ki), the above lemma reduces to Lemma ^ 

Proof Again we will rely on the induction argument. As it was already mentioned, the general soliton 
matrix T(k) corresponding to several distinct zeros can be represented as a product (|45|) of individual 
soliton matrices (J27|) for each zero. For clarity reason and simplicity of the presentation we will give 
detailed calculations for the simplest case of just one product in (|45jl . Then we will show how to 
generalize the calculations. Consider soliton matrix T(k) for two pairs of distinct higher-order zeros 
(k±,ki) and (&2, ^2)- We have T(k) = Ti(k)T2(k) and 

T(k)=(l + -^ + ...+ , A _\ Wl + -g^ + ,,, + - ggi- V (49) 



k-ki (k-ki) n ^J\ k-k 2 (k - k 



"2, 



11 2 



Here n 3 - (j = 1, 2) is the number of simple matrices in the product representation (J27j) for Lj. Due to 
Lemma ^ the coefficients Aj and Bj are given by formulae similar to (|36p: 



r!" .<'>+l-j 



i/=l Z=l 



S »=E E i^WU-J- (5i) 

i/=i i=i 

On the other hand, by expanding formula (f4H|) into the partial fractions we get 



r(jfc) = 1 + + ...+ - + -i- + . . . + - ^ — . (52) 

fc-Jfel (A; - /ci) ni A;-A; 2 {k - k 2 ) n * 

Consider first the coefficients Aj. Multiplication by (k — ki) ni of both formulae ()49j) and ()52|1 and 
taking derivatives at k — k\ using the Leibniz rule gives 



In similar way we get 

*.- = £^&)^r («) 

3=0 v J; 

Now substituting Eqs. (fBTljl and (fBTj) into ((HSj) and (JH4"|) and defining new vectors 

1 m ~ 1 1 d-T 

<Pm I = JZ^-il^TfcT^ 1 )' m = V--,^. (55) 

j=0 J ' 

and 

m— 1 1 

~(u,2) ^ 1 d^Ti - 



\Qr. 



> = £^&)l&3>. m =!,...,«?>, (56) 



3=0 J 
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we find that 



i/=i i=i 



+2-j -V 



sir'+2-j-r' 

i/=l Z=l 

which give precisely the needed representation (J46j) . Note from definitions (J55j) and (J56|) that 



(57) 



(58) 



and 



(Pi 



p(^2) 



(P r (D I 
' i 



,~(^2) ' 
l?„(2) ) 



Ti(fc 2 



i# 2) ) 



r 2 (fci 



iy (2) / 



Due to lemma HI vectors {(Pi |,---,(p^i^|} and ), • • • , |y ( 2 ) )} are linearly independent 



respectively. 



are non-degenerate. Thus new vectors 



In addition, matrices T 1 (k 2 ) and r 2 (fci) 
{(pi ' ; |, • • • , (p\i) 1} and {l^i ' ),..., I? (2) )} are linearly independent respectively as well. This com- 
pletes the proof of Lemma [2] for two pairs of higher-order zeros. 

It is easy to see that the above procedure of redefining the vectors in the blocks corresponding to 

different zeros will also work in the general case, when Ti(k) is replaced by the product Ti(k) r n (A;), 

and T 2 (k) replaced by T n+ i(k). In this case, the sum over all distinct poles will be present in the left 
Q-bracket in formula (J49j) . and consequently there will be more terms in formula (J52j) . Formula (J53|) 
will be valid for coefficients A of each zero, and formula ()54j) remains valid as well. Thus by defining 

vectors (jp^f \ by formula for each zero kj (1 < j < n), and defining vectors \q m + ) by formula 
()56|) for zero k n+ i, we can show that the matrix T(k) consisting of n + 1 products of Tj(k) can be put 
in the required form (j4"Hjl . This induction argument then completes the proof of Lemma EJ Q.E.D. 



The notations in the representation (|46|) for soliton matrices with several zeros are getting com- 
plicated. To facilitate the presentations of results in the remainder of this paper, let us reformulate 
the representation ()46j) . For this purpose, we define r\ 



+ ... + r 



(Nz) 



where r^'s are as given 



in Lemma 121 Then we replace the double summations in Eq. (J46|) with single ones, 



r(*) = / + J2 s v , r- 1 (k) = i + J2 s - 



(59) 



v=l 



v=l 



Inside these single summations, the first terms are blocks of type (J47j) for the first pair of zeros 
(k\, ki), the next r| 2 ^ terms are blocks of type (jUj) for the second pair of zeros (kz, k 2 ), and so on. 
Block matrices S u and S u can be written as 



4 = EE 



1 = 1 j=l V V > 



-/ 



'(pfv 



(60a) 
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' I,,"" \ /,,(") I 



c 1-PiVi-jX* 



1=1 j=l y v > 



+1-1 



V<rf°iy 

where matrices D v (k) and D v {k) are triangular Toeplitz matrices with poles: 



(60b) 



/ 



i 



1 





1 



\ 



/ 



, A,(fc) 



i i i 



1 



1 



1 







1 



1 







1 



(61) 



(k-K v ) J 



Here 



k 



3-1 

^ .(0 
1 



3 



if l + }_^r\ l > <v<J2 r ? 
i=i i=i 



;i < 3 < N z ] 



(62) 



In other words, k v = k\ for 1 < v < , = &2 for rj 1 ^ + 1 < z/ < + r\ , etc. In addition, 
{s„, 1 + 52j=i — v — Yl{=i r i } i s the block sequence of the j-th pair of zeros (kj, kj). This new 
representation (jSUjl is equivalent to (JIB]), but it proves to be helpful in the calculations below. 

We note that the economical way of block numeration used in the representation (}59|) reflects the 
important property of the solitons matrices: the soliton matrices preserve their form if some of the 
zeros coalesce (or, vise versa, a zero splits itself into two or more zeros). The only thing that does 
change is the association of a particular z/-block to the pair of zeros. 

The representation [or (|3B))] is but the first step towards the necessary formulae for the soliton 
matrices. Indeed, there are twice as many vectors in the expressions (15 9 j) for T(k) and r _1 (A;) as 
compared to the total number of vectors in the constituent matrices in the product of representations 
of the type (}2*T|) for each pair of zeros. As the result, only half of the vector parameters, say \p^) 
and (p^ |, are free. To derive the formulae for the rest of the vector parameters in (|59J) we can use 
the identity r(fc)r _1 (/c) = Y~ l {k)Y{k) = I. First of all, let us give the equations for the free vectors 
themselves. 



Lemma 3 The vectors \p^), ■ ■ ■ , \pfj) and {p( 
( E^ -( En|) satisfy the following linear systems of equations: 

( 



,---ApV\ 



TJk v ) 



( \pV) ^ 



\ \pt ] ) J 



o, r„(fc) = 



1 _d_ 

1! dfc 



from each v-th block in the representation 

\ 

(63) 





r 



o 



o 



1 d 
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Remark Note that the matrices T u (k) and T u 1 {k) have block-triangular Toeplitz forms, i.e., they 
have the same (matrix) element along each diagonal. 

Proof The derivation of the systems (fM|) - (j64j) exactly reproduces the analogous derivation in Ref. 
for the case of elementary zeros (as the equations for the z/-th block resemble analogous equations 
for a single block corresponding to a pair of elementary zeros). For instance, the system (}63|) is 
derived by considering the poles of T{k)Y~ 1 {k) at k = k, u , starting from the highest pole and using 
the representation (JoT?j) - (JoT?|) for T^ 1 (k). The details are trivial and will not be reproduced here. Note 
that there may be several sets of vectors (from different //-blocks of the same pair of zeros) which 
satisfy similar equations if the geometric multiplicity of this pair of zeros is higher than 1. Q.E.D. 

Now let us express the \q)- and (g|-vectors in the expressions (|59|) - (j6lH) for T(k) and r _1 (A;) through 
the \p)- and (p|-vectors. This will lead to the needed representation of the soliton matrices given 
through the \p)- and (p|-vectors only. It is convenient to formulate the result in the following lemma. 

Lemma 4 The general soliton matrices for several pairs of normal zeros {(ki, k\), . . . , (k^z, k^z)} 
are given by the following formulae: 




T(k) = I-(\p?) 



(65a) 



(Pi 




<pg?l 



T-\k) = I - (\p?) 



...,b«),...,br ) ),...,|pW)MA:)/C 



-i 



(65b) 



(ft 
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where s u and n are the same as in Lemma\^ The matrices T>(k) and T>(k) are block- diagonal: 



1 D 1 (k) 



V(k) = 



\ 







\ 



D ri (k) ) 



( 



Vik) 



D 1 (k) 



\ 







\ 



D ri (k) ) 



(66) 



where the triangular Toeplitz matrices D u (k) and D u (k) are defined in formulae \61)) . The matrices 
K and K have the following block matrix representation: 



K, = 



with the matrices K v>tt and K^'^ being given as 



K 



( . . . #( 1 ' r \ 



(67) 



K 



EE 

3=0 1=0 



1=0 3=0 



Here {H_] , 1; . . . , Hg"!_i} is the basis for the space of s v x s v - dimensional Toeplitz matrices, i.e., 

(Hj"^) at p = $ a ,i3-j- The nonzero elements of matrices Qi"'^ and Q\ u '^ are defined as the inner 
products between the p-vectors from the blocks with indices v and \x: 

( \ 



r(") 



Hi — 



\¥2\/ 



(0,...,0,|^),...,bi^)), Q\ 







Remark 1 In the case of a single pair of zeros (ki,ki), we simply replace (k^) and k v (k v ) in 
formula (|6T|) by k\ {k\). 

Remark 2 In the case of the involution (jJJ property, we have the obvious relations: 



(|rfV-. >£?>). (69) 



Proof We only need to prove that the \q) and (g| vectors in soliton matrices (j59jl -(|60 jl are related to 
the \p) and {p\ vectors by 



(|g«), . . . , |^>, • • • , |g! n) ))/C = -(\p\ l> ), |pW>, • • • , IpD, - - - , K J )), (70) 



(i)\ 



>0\ 



,(n)\ 
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and 



( (p?\ \ 



K. 



<<£> 
fa 



(n)| 



(M n) l 



(71) 



V^ ri) l/ 



where matrices /C and /C are as given in Eq. (|H7j) . We will give the proof only for Eq. (fT0|) , as the 
proof for (j71|) is similar. Note that in the case of involution (J3J), Eq. (j71|) is equivalent to ()70|) by 
taking the Hermitian. 

To prove Eq. (J7TJ), we consider the corresponding expression (loTIJl - lJoTH) for T(k): 



Y{k) = I + Y J i\t]),---Mi ) ))DM 



(72) 



We need to determine the |g)-vectors from Eq. (JHSI)- Note that the /-th row in the //-system (IFk?|) can 
be written as 

, ,i dr i d^r 

rM ' n dT ( ^ } ' • • • ' (z-i)!dfc<- iKt) 







(73) 



for each 1 < \i < r±. When the T(fc) expression (fT2~j) is substituted into the above equation, we get 



ri 



1 (#1 ^ 



Ed^ ) >.---.i9S" , ))i 



1! dfc 



IpS) 



+ ••• + - 7TT 7TTTI (*W 



(/-!)! dk l ~ 



\ 



\p { r ] ) 



> = -irf°>. 



(74) 



The derivatives of D v (k^) can be easily computed: 



d'D ^ (-lyq+Q i #17 



(75) 
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Now it is straightforward to verify that all equations of the type (j74j) can be united in a single matrix 
equation (|70j) by padding some columns in the summations of (|73jl by zeros, precisely as it is done 
in the definition of Q . As a result we arrive at the relation (|70j) between \q) and |p) vectors, 
where the matrix K, is precisely as defined in Lemma 0J Q.E.D. 



D. Two special cases 

The soliton matrices derived above reproduce all previous results as special cases. Previous results 

'eometric 
22|. In the 



were obtained in two special cases: several pairs of Riemann-Hilbert zeros with equa 
and algebraic multiplicities (lif . and a single pair of elementary Riemann-Hilbert zeros 
first case, suppose the geometric and algebraic multiplicities of n pairs of Riemann-Hilbert zeros 
{(kj, kj\ 1 < j < n} are {r^\l < j < n\ respectively. Then the soliton matrices have been given 
before [13| (see also appendix B in Ref. |35|) as: 

n P W r W) | Wwn-l\ /st(0| n r« pO) , (l)wp-U /-Mi 

r = J -EEE - ' r " = / + EEE - Jr ' (™> 

j j=l m=l 1=1 J i ,j=l m=l 1=1 ^ 

where vectors {|i>j ), 1 < / < r^} and 1 < Z < r^} are in the kernels of T(kj) and F~ 1 (kj) 

respectively: 



r(%)k (/) ) = o, ^|r- 1 (* i ) = o, / = l,...,rW (77) 



and 



Moreover, 



im,jl 



det r = J~{ 



/_(m)| (0\ 



A? 



(78) 



,0) 



3 =1 \ 3 

The above special soliton matrices can be easily retrieved from the general soliton matrices (|65|) -(|69| ) 
of lemmaEJ Indeed, in this special case, the block sequence of a pair of zeros (kj, kj) is consecutive 
l's. Thus s v = 1 for all z/s. Consequently, matrices D v and D u in Eq. ()66|) have dimension 1. In 
addition, matrices K^ v ^ and in Eq. 1)68)1 also have dimension 1, and the summations in their 

definitions can be dropped since I = and j = there. Hence, we get 

j^y,n) = r K (»,v)y = (P i i ) \Pi l) ) ^ 

see (JEHJ) - Relating |p)-vectors {|p^ } >, 1 + JXjrW < v < ELi r(0 } to {I^M < 1 < rU) } and 
{(p[ \, 1 + Yli=i r ^ — v — Y^i=i ^° i( v j I) 1 — ^ — for each j = 1, . . . , n, and recalling the 
definition of /t's, we readily find that our general representation reduces to (JTHJ). We note 
by passing that the soliton matrices (|7nj) - (|7Hj) cover the case of simple zeros, where there is just one 
vector in each kernel in ()77j1 . 
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Our second example is a single pair of elementary higher-order zeros. A higher-order zero is 
called elementary if its geometric multiplicity is 1 [22]. This case has been extensively studied in the 

literature before (see Refs. EE0EEE3) 

for different integrable PDEs. The soliton matrices having 
similar representation as ([65 )1 - ([69 )1 for this case were derived in our previous publication j^. The 
only difference between that paper's representation and the present one ([65 p -(|69| ) is the definition 
of the matrices /C and /C. However, in this special case, these matrices have just one block each, 
i.e., K^' 1 ) and K 1,1 , since there is just one z/-block in the soliton matrices. By comparison of both 
definitions one can easily establish their equivalence. 



E. Invariance properties of soliton matrices 

In this subsection, we discuss the invariance properties of soliton matrices. When the soliton matrix 
is in the product representation (}2Tj) for a single pair of zeros, the invariance property means that one 
can choose any r 1 linearly independent vectors in the kernels of r(ki) and r _1 (A;i), or more generally, 
one can choose any r\ (1 < I < n) linearly independent vectors in the kernels of (IAj -1 ■ ■ ■ xT-i)(^) 
and (xi-i ■ ■ ■ XiT ^(^l), and the soliton matrix remains invariant. In other words, given the soliton 
matrix T(k) for a fixed set of 77 linearly independent vectors \vu) (1 < i < 77) in the kernels of 
(IA^ 1 . . -xT-i)^) and another fixed set of 77 linearly independent vectors (vu\ (1 < i < 77) in the 
kernels of (xi-i ■ ■ ■ Xi^ 1 )^!), new sets of vectors 

[\vu), \v 2 i), . . . , \v n ,i)] = [K), \v2i), • • • , \v n ,i)\ B, (79) 

and 



( v n\ 




( v u\ 


(V2l\ 


= B 


(V2l\ 






(v n ,i\ 



(80) 



where B and B are arbitrary /c-independent non-degenerate r\ x r\ matrices, give the same soliton 
matrix T(k). This invariance property is obvious from definitions (|25|) for projector matrices. Note 
that the invariance transformations (fT9"j) - (J%Uj) are the most general automorphisms of the respective 
kernels (null spaces) (r^ 1 . . . x7-i)(ki) an d (Xi-i ■ ■ ■ Xi^ r )(ki)- 

Now let us determine the total number A/f ree of free complex parameters characterizing the higher- 
order soliton solution. For a single pair of the higher-order zeros (ki, k\) in the case with no involution, 
it is given by the total number A/tot (= 2A^" =1 r^ + 2) of all complex constants in all the linearly 
independent vectors in the above null spaces and the pair of zeros (k\, ki), minus the total number 
A/inv (= 2 Ym=i r f) °f ^ ne f ree parameters in the invariance matrices (|79*|) -(|5n |) . Thus, in the case with 
no involution, we have 

n n 

Mree = Mot - Ai nv = 2 A ^ + 2 - 2 ^ t] . (81) 

1=1 1=1 

Note that the total number of \v) or (v\ vectors in the product representation ([27)1 . given by the sum 
Y^i=i r u is equal to the algebraic order of the pair of zeros (ki, k{). In the case of the involution (j3J), 
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the number A/f ree is reduced by half. When the soliton matrices have several pairs of zeros as in the 
product representation ([45 jh the invariance property is similar, and the total number of free soliton 
parameters is given by the sum of the r.h.s of formula (|81j) for all distinct pairs of zeros. 

By analogy, the invariance properties for the summation representation (J55"j) of the soliton matrices 
are defined as preserving the form of the soliton matrices as well as the equations defining the \p)- and 
(p|-vectors (|63| ) -([64 p . The equations defining the transformations between different sets of p- vectors 
of the same invariance class must be linear, since all the sets of p-vectors in the invariance class 
satisfy equations (|S3 )) -(|B4~ }) for a fixed soliton matrix - i.e. the invariance transformations are a subset 
of transformations between solutions to a set of linear equations. Thus the most general form of the 
invariance is given by two linear transformations — one for |p)-vectors and one for (p|-vectors: 



lfff>, ■ ■ • ■ IpS"'), ■ ■ ■ ■ |p£'» = • • • , Ip2'), • • • . \pt' ) ), \P { I r \'))B, (82) 



and 



(pi 



In), 



\ /-( ri >\ I 



I (p?\ \ 



(p£>\ 



B 



(n)i 



(83) 



Different from the product representation of the soliton matrices, the transformation matrices B and 
B in Eqs. f!82j) and (|H3*jl can not be arbitrary in order to keep the soliton matrices and equations 
invariant. Let us call such matrices B and B which keep the soliton matrices (ffiKj) invariant 
as invariance matrices. The forms of invariance matrices can be determined most easily by considering 
the invariance of equations (j53 |) -(p4" |) . 

Recall from Lemma El that all \p) vectors in the soliton matrix (jUB*]) satisfy the equation 



ri(«o 



(84) 



( \P?) \ 
\P$) 
lrf ri) ) 

V [pfe } > J 

where T v (k v ) is the lower-triangular Toeplitz matrix defined in Eq. ()63|). The matrix B is an 
invariance matrix if and only if the above equation is still satisfied when the \p) vectors in Eq. ()84j) 
are replaced by the transformed vectors \p) in Eq. (|82j) . and the resulting matrices K, and /C are 



r ri (/t ri ) j 
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non-degenerate [see Eq. (|65|)]. Note that the transformation (|82|) can be rewritten in the following 
form: 



( m \ 

iM ri) ) 



B 



T 



( \P?) \ 



lM n) ) 



(85) 



where the superscript "T" represents the transpose of a matrix. Since the original \p) vectors can be 
chosen arbitrarily (the matrix Tb is determined subsequently from these \p) vectors as well as the 
(p\ vectors), in order for the above \p) vectors (jHSJ) to satisfy Eq. (|51j) as well, the necessary and 
sufficient condition is that Tb and B T commute, i.e., 



B 1 



B 1 



(86) 



and B is non-degenerate. The requirement for the non-degeneracy of B is needed in order for the 
resulting matrices /C and /C to be non-degenerate [see Eq. (|55jl]. Similarly, we can show that the 
matrix B in Eq. (|H3|) is an invariance matrix if and only if Tb and B commute, 



T n B 



b t -r B , 



if) 



and B is non-degenerate. Here the block-diagonal matrix Tb is 

\ 



\ 







(88) 



and upper-triangular Toeplitz matrices T v (k u ) have been defined in Eq. (jMJ). Note that matrices T B 
and Tb have exactly the same forms as T>(k) and T>(k) respectively. Thus invariance matrices B T 
and B commute with T>(k) and T>(k) as well: 



V(k) -B 1 =B L -V(k), 



V(k)-B T = B T -V(k). 



(89) 



In addition, since T) T has the same form as V, invariance matrices B and B also commute with T> 
and D: 

B ■ V{k) = V{k) ■ B, B ■ V(k) — V(k) ■ B. (90) 

The forms of these invariance matrices are easy to determine. First of all, the commutability 
relations (J9U|) demand that the invariance matrix B has a block-diagonal form with each block cor- 
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responding to a pair of zeros: 



/ 



B 



\ 



Bo 



\ 



B 



(91) 



Nz J 



Here B n is a square matrix associated with the n-th pair of zeros (k n , k n ). The form of each matrix 
B n is readily found to be 



Dn ■ ■ ■ D n 



B r 



R (ri n) ,l) R (ri n) ,rW) 
-t>n • • • -Dn 



(92) 



where Bn is a sl™^ x matrix of the following type: 











6l 6 2 ••• 6,(n)_i bin) 



6 X fe 2 



6 .< B) -i 



h b 2 
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U>fi, 



(93a) 



Ci c 2 ... 

ci c 2 



V 



... 



cm 
c 2 

Cl 






U<fl, 



(93b) 



/ 



sS n) > 4 n) > • • • > is the block sequence of zeros (k n ,k n ) as in Lemma El (see Definition EJ), and 

6j,Cj are arbitrary complex constants which are generally different in different submatrices Bn'^. 

The invar iance matrix B has the form of B T . 

The above forms of the invariance matrices B n and B n follow immediately from the 

following argument. Consider, for instance, the matrix B n . The commutability relation with the part 

of the matrix T>(k) corresponding to the ra-th pair of zeros, i.e., T>^ n \k) = diag[.D^(/c), D^l^k)] 

r i 

where matrices D„ (k) are given by Eq. (JIEJ), produces the following set of independent matrix 
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equations 

Df\k)B^=B^D^\k), 1/^ = 1,...,^. (94) 

For v = /i, the above equations are equivalent to the commutability conditions for the single elemen- 
tary higher-order zero considered in Ref . (2^ , thus the form (J93|) for the diagonal blocks Bn^ follows 
accordingly. Consider now the case when v > \i (the other case can be considered similarly). We 
have then s£ < and the square matrix \k) contains the matrix Dy(k) in its lower right 
corner [consult the definition ([48)1] . It is easy to conclude, first of all, that the first fi — v columns of 
the matrix Bn are identically zero, otherwise on the r.h.s. of equation (j94j) we would have higher 
powers of (k — /c„) _1 than the highest power of this quantity on the l.h.s.. Then if we denote the 
non-zero part of Bn'^ as Bn \ the condition (}9*4*j) becomes 

DS>\k)BM = M*DM(k), 

which is equivalent to the one considered above in the case of \i = v. Thus the form (J93b ) for the 
off-diagonal blocks of the invariance matrix Bn follows as well. Q.E.D. 

From the above explicit expressions (|91j) -()93 j) for invariance matrices in the summation represen- 
tation (JOS}; it is easy to see that the total number A/l nv of free complex constants in these invariance 
matrices coincides with that in the product representation pTj) and (}45|) [see Eq. (J8lf ]. Indeed con- 
sider for simplicity just a single pair of zeros. In the case with no involution (0J), the total number 
A/lnv of free complex constants in the invariance matrices (}9"Tj) -(f9"3 j) is 

n n 
Af inv = 2 Oi - 2z/ + 1) s ri ^ +1 = 2 J2{2fj, - l)s, 



(r n r n -i n \ n 

n£(2Ai-l) + (n-l) ^ (2/x - 1) + ... + ]T ( 2 M " 1) = 2 
fi=l /j,=r n +l fi=r 2 +l / 1=1 



(95) 



which is exactly the same as that in Eq. (|%T|) for A/lnv ■ Here we have used the fact that the numbers 
of blocks with sizes [1, 2, 3, n] are given by the differences of the ranks [ri — r 2 , r 2 — r 3 , r n _i — 
r n, r n] (see end of Sec. IIIIB|) . This result is not surprising since the invariance properties of the 
soliton matrices in the summation representation originate from the invariance properties in the 
product representation, that is why the respective invariance matrices have the same total number 
of free parameters. Consequently, the total number of free complex parameters in the summation 
representation (|65j) is the same as in the product representation, in the case with no involution for a 
single pair of zeros it is given by the same Eq. ([81)1 . 

Invariance matrices have many important properties. These include (i) the identity matrix / is 
an invariance matrix; (ii) if B is an invariance matrix, so is cB, where c is any non-zero complex 
constant; (iii) if B is an invariance matrix, so is -B -1 ; (iv) if B\ and B2 are two invariance matrices, 
so are B\ ± B2 and B\ ■ B2. In the former case, B\ ± B2 should be non-degenerate. 

Lastly, we note that if matrices B and B satisfy the commutability relations ([90)1 . the transforma- 
tions (|82j) and (JHHj) indeed keep the soliton matrices invariant. The proof uses the fact that under 
the transformation ()82j) where B is an invariance matrix (the (p| vectors are held fixed), matrices K, 
and /C are transformed to _ 

jt = KB, K = KB (96) 
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respectively. Similarly, under the transformation (|83|) while keeping the \p) vectors fixed, matrices /C 
and K are transformed to _ 

jC = BJC, K = BK. (97) 

For a single pair of elementary zeros, these facts have been proved in j^. The proof for the present 
general case is given below. Since the proofs for Eqs. (}9l)j) and (}9*Tj) are similar, we only consider Eq. 



To prove the transformation (J96|) . we need to recall how matrices K, and /C are obtained. The 
matrix JC is derived from Eq. (|84j) . Comparing this equation with (|7()jl . we find that 



I) 



\p$) 
\p[ ri) ) 



— T 

K 



( m \ 



tit?) 



V l^ ri) ) / 



Now under the transformation (|82)l. i.e., (|85j) . and recalling that B T and Tb — / commute, we readily 

— T — T — — 

find that /C = B T K , thus /C = /C-B. As about the matrix /C, it is derived from the equation 

(^|,... ) <pg)|,... ) ^|,... ) ^)|)r B = o, 

where T B is given by Eqs. (Jfi3|) and (J%%|) . Recall that T~ 1 (k) is given by Eq. (jKHjl . i.e., 

/ fog 3 ! \ 



r- 1 (fc) = / + (|ri 1) ),...,bi 1 i )),. 



|pS n) ),- 



(9i 



(n). 



(V 1 



\(i{ ri) \J 

Under the transformation (|82jl . noting that .B and T> commute [see Eq. (|9()j)]. we readily find that 
K = KB. Thus flllj) holds. 

Because of Eq. (|96p and the commutability relation (|90|). we see that soliton matrices r(/c) and 
T _1 (A;) in Eq. ()65|) indeed remain invariant under the transformation (|82|) . Analogously, these soliton 
matrices are also invariant under the transformation (J%3|) if matrix is an invariance matrix. In 
the case of involution (jlj), transformations (JH2|) and (JH3|) need to be performed simultaneously since 
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\p) and {p\ vectors are related by the Hermitian operation. Under these combined transformations, 
matrix /C transforms to /C = BICB, thus soliton matrices ()65|) remain invariant as well. 

The invariance matrices can be used to reduce the free parameters in soliton solutions to a minimum 
as we have done above [see Eq. (|%T])]. It is also needed to classify the general evolution of soliton 
matrices (see the next subsection). 



F. Spatial and temporal evolutions of soliton matrices 

Finally, we derive the (x, ^-dependence of the vector parameters which enter the soliton matrix 
(foT)|). The idea is similar to that in the derivation of equations (J2U|) in section |TT1 Our starting point 
is the fact that the soliton matrix T(k,x,t) satisfies equations ©-© with potentials U(k,x,t) and 
V(k,x,t): 



d x T(k,x,t) = T(k,x,t)A(k) + U(k,x,t)T(k,x,t), 
d t T{k,x,t) = T(k,x,t)Q{k) + V(k,x,t)F(k,x,t). 



(98a) 
(98b) 



First we need to find the equations for the triangular block- Toeplitz matrices T v and T u . To this goal 
one needs to differentiate equations (}9"%j) with respect to k up to the (s u — l)-th order. It is easy to 
check that the equations for T u have the same form as equations 



d x T u (k, x, t) = T v {k,x,t)h v {k) + U u (k,x, t)T v (k,x, t), 
d t T u (k,x,t) = T v {k,x,t)n v {k)+V v {k,x,t)T v {k,x,t). 



(99a) 
(99b) 



Here A u , Q u , XJ U , and V„ are lower-triangular block- Toeplitz matrices: 
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Indeed, this is due to the fact that the matrix multiplication in (J99)) exactly reproduces the Leibniz 
rule for higher-order derivatives of a product. Similarly, using the equations for T -1 , 



d x T-\k,x,t) 
d t T-\k,x,t) 



-A(Jfe)r -1 (A;, x, t) - T-\k, x, t)U{k, x, t), 

-nC^r-^ife, x, t) - r- x (A;, x, t)v(k, x, t), 



(102a) 
(102b) 



28 



one finds that 

d x T u (k,x,t) = -A v (k)T v (k,x,t) -T v (k,x,t)V v (k,x,t), 
d t T v (k,x,t) = -fl„(k)r v (k,x,t) -T v (k,x,t)V v (k,x,t), 

where A„, fi„, XJ U , and *V U are upper-triangular block- Toeplitz matrices: 
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(103b) 
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To obtain the (x, ^-dependence of the p- vectors, we differentiate equations (|U3*j) and Utilizing 
Eqs. (HUTO and (fTU^ . we find that 



[9a, + A„(«„)] 



(106) 



and 



T v {k v ) 



[d t + n v { Ku )} 



> = o. 



(107) 



Due to the invariance properties (see explanations below), we can set the quantities inside the curly 
brackets of Eqs. ()106|) and ()107|) to be zero without any loss of generality: 



[d x + A„(«„)] 



o, [dt + n u {n v )} 



(108) 



The reason for it is the uniqueness of solution to the Riemann-Hilbert problem for the given spectral 
data. Thus, the (x, Revolution of \p) vectors is 



( \vf) \ 



1 \vth ^ 



exp{-A iy (^)x - Vt v {K v )t\ 
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(109) 
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By similar arguments, the (x, ^-evolution of 



vectors is 



(110) 



Here the subscript "0" is used to denote constant vectors. The exponential functions in the above 
two equations can be readily determined. Indeed, by using the property that the operation of raising 
to the exponent of a diagonal matrix (such as A(k)x + Q(k)t here) commutes with the construction 
of the block- Toeplitz matrix (see appendix in Ref. we find that 
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where E{k) = exp {— A(k)x — Q(k)t}. After the spatial and temporal evolutions of vectors \p) and (p\ 
have been given from Eqs. (jl09|) to (jlll|) . the soliton matrices f!65|) are then obtained. Eventually, the 
soliton solutions are derived from Eq. ((SJ) by taking the limit k — > oo. For the three- wave interaction 
model, soliton solutions are given by Eqs. ()12|) and (fTSj) . The corresponding eigenfunctions of the 
N- dimensional Zakharov-Shabat spectral problem with those soliton (reflection-less) potentials are 
simply the column vectors of the soliton matrices T(k) and T~ 1 (k) in (J65|) by taking k to be zeros 
(kj, kj) (which give discrete eigenfunctions) and with k lying on the real axis (which give continuous 
eigenfunctions). 

Lastly, we explain why other \p) solutions to Eqs. (jK)fij) and ()107j) give the same soliton matrices 
as those from Eq. (jl08)l . Notice that equations (jlOfij) for all v blocks can be written in the following 
compact form: 

/ \P?) \ 
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According to the invariance properties in the subsection IIIIE1 any two vectors in the kernel of matrix 
Tb are linearly dependent. Thus the most general \p) solutions to Eq. f)106j) are such that 



lM n) ) 
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|jfi } > 
|Pi n) ) 
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where B is an invariance matrix which depends on x and t in general [see Eq. To show that 

these \p) vectors give the same soliton matrices (jp^|) as the \p) vectors from Eq. ()108|) . we define a 
matrix function G(x,t) which satisfies the following differential equation and initial condition: 



8 x G{x,t) = B T {x,t)G{x,t), G\ x=0 



I. 



Because the matrix B here is an invariance matrix and G(x = 0) = I, obviously the function G(x,t) 
is an invariance matrix as well (note that G is always non- degenerate from its construction). In 
addition, G~ l is also an invariance matrix. Now for any solution \p) of Eq. ()112|) . we define new 
vectors \p) as 



( \P?) \ 
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G- 1 



\ \p£?) ) \ J 

Then these \p) vectors satisfy the first equation in p08|) . This can be checked directly by substituting 
the above equation into (|1U8|) and noting that matrices G and fie commute by virtue of Eq. (|HU|) 
and the fact that matrices Qb and Tb have identical forms. Since G~ x is an invariance matrix, \p) 
and \p) vectors as related above naturally give the same soliton matrices (|65j). Thus there is no any 
loss of generality by picking solutions \p) of Eq. (jl(J6j) such that the first equation in f)108j) holds. By 
the same argument, there is no loss of generality by picking solutions \p) of Eq. (J1U7)) such that the 
second equation in (jl()8|) holds. 



IV. APPLICATIONS TO THE THREE- WAVE INTERACTION SYSTEM 

To illustrate the above general results, we apply them to the three-wave interaction model 
and display various higher-order soliton solutions. In this case, the involution property (jlj) holds, thus 
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all zeros are normal and appear in complex conjugate pairs. The soliton matrix T(k) is given by Eq. 
(|55k). where {p\ = \p)\ and the (x, ^-evolution of \p) vectors is given by Eqs. ()109|) and (|lllh ). The 
general higher-order soliton solutions of the three- wave system are then given by Eq. ()13j) . where 



$w = r« = -(b«),...,b«),...,br i} ) 



(ri)\ 



\pt ] ))K 
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\(p { ^\J 

and matrix K is given in Eq. (|67|). In all our solutions, we fix the parameters in the dispersion laws 
(ITT1) as (ai,a 2 ,a 3 ) = (1,0.5,-0.5) and (61,62,63) = (1,1-5,0.5). 



(114) 



A. Soliton solutions for a single pair of non-elementary zeros 

First, we derive soliton solutions corresponding to a single pair of non-elementary zeros. In partic- 
ular, we consider the rank sequence {1, 2} of a pair of zeros (k±, k\). In this case, r\ — 2 and r 2 = 1. 
Using formula (JHTj) (for the case of involution) we get the number of free complex parameters in the 
soliton solution: 

TVfrce = 3(2 + 1) + 1 - (4 + 1) = 10 - 5 = 5. 

There are three \p) vectors, \p x ), |?4 ) and \p^) in Eq. ()114|) . When k\ and the initial values 
[boi/j I.P02 )> \Pm)] °f t nese vectors are provided, the soliton solutions (fT3"j) will then be completely 
determined. 

In the present case, the block sequence reads {s\, s 2 } = {2, 1}, and the invariance matrix B for 
this case can be readily obtained from the general formula (fHTj) as 
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which indeed has five free complex parameters [see Eq. (|95J) ]. The invariance matrix B is just the 
Hermitian of the B matrix. 

To display these soliton solutions, we choose k% = 1 + bo^) = [ — 1, 1 — boi ) = [1, 0.5, — l] r . 
When \p$) = [l,l + i,0.5] T (the generic case), the solutions are plotted in the top row of Fig. 1. 
In two non-generic cases (where some elements of the \p) vectors vanish), \p$) = [0, 1 + z,0.5] T and 
boi ) = [1; 0; 0.5] T , the solutions are plotted in the second and third rows of Fig. 1 respectively. We see 
that in the generic case, three sech waves in the three components interact and then separate into the 
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same sech waves with their positions shifted. In other words, this is a wi(sech) + w 2 (sech) + w 3 (sech) — > 
tii(sech) + w 2 (sech) + w 3 (sech) process. What happens is that the initial pumping (w 3 ) wave breaks 
up into two sech waves in the other two components (u% and u 2 ), while simultaneously the two 
initial U\ and u 2 waves combine into a pumping sech wave. Thus this process is a combination of 
two sub-processes: w 3 — ► U\ + u 2 and U\ + u 2 — > w 3 . This phenomenon seems related to the rank 
sequence {1, 2} of the present solitons and the fact that, the rank sequence {1} itself describes 
the breakup of a pumping sech wave into two non-pumping sech waves, while the rank sequence 
{2} itself describes the reserve process. In the non-generic cases, these solutions can describe the 
«i (sech) + u 2 (second-order) — > w 2 (sech) + w 3 (sech) process, the Wi(sech) + w 2 (sech) + -u 3 (sech) — > 
m 3 (second-order) process (see Fig. 1, second and third rows), and many others. In the solutions of 
Fig. 1, the dj and bj parameters are such that u 2 < 1*3 < u\. If U\ < u 3 < u 2 , the processes will be 
exactly the opposite (see 2^1 ). Thus our solutions can describe the opposite processes of Fig. 1 as 
well. 



B. Soliton solutions for two pairs of simple zeros 

Here we derive soliton solutions corresponding to two pairs of simple zeros in the three-wave system 
(I14|) . Some solutions belonging to this category have been presented in 0, |2t|. But we will show 
that those solutions are only special (non-generic) solutions for two pairs of simple zeros. Below, the 
more general solutions for this case will be presented. 

In this case, = = 1. By using formula (JHTj) for the case of involution (j3J) and with two 
pairs of zeros, we readily obtain that the number of free complex parameters in the solution is 6: 

jVfree = 2(3 X 1 + 1 - 1) = 6. 

Indeed, there are two \p) vectors in Eq. (jll4J) . Together with the two zeros k\ and k 2 , there are 8 
complex parameters in the soliton solutions. However, the 2x2 invariance matrix B in this case is 
diagonal and has two free (diagonal) complex parameters. 

Three solutions, with k\ = 1 + i, k 2 = — 1 + 0.5z and three different sets of |j>oi ) an d bo?) vectors, 
are displayed in Fig. 2. In the generic case where \p^i) = [1,1 -M,0.5] T and |pqi ) = [1, 0.5, — l] r 
(see top row of Fig. 2), the solution describes the breakup of a higher-order pumping (w 3 ) wave 
into two higher-order U\ and u 2 waves. This is analogous to solutions for a single pair of elementary 
zeros with algebraic multiplicity 2 (see J22j). In the non-generic case where \p$) = [0,1 +z,0.5] T 
and |j»oi ) = [1)0.5,— 1] T (second row in Fig. 2), the present solutions can describe the w 2 (sech) + 
w 3 (sech) — > ui(sech) + n 2 (second-order) process. This process has been seen in jij] for elementary 
zeros as well. More interestingly, in the non-generic case when P01 [1] = P01 [3] = 0' these solutions 
describe the elastic interaction of a sech u\ wave with a sech u 2 wave (see bottom row of Fig. 2). 
These are precisely the soliton solutions presented in 2(| 2?| . We see that these solutions are simply 
non-generic solutions for two pairs of simple zeros. 
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C. Soliton solutions for two pairs of higher-order zeros 



Lastly, we consider two pairs of zeros, one simple and the other one elementary with the algebraic 
multiplicity 2. Let us say k\ is the elementary zero, and k 2 (7^ k\) is the simple zero. Then the rank 
sequence for k\ is {1, 1}, and the rank sequence for k 2 is {1}. Thus, 
By formula (JHTj) we have 



1, = 1, and r^' = 1. 



(2) 



free 



3(l + l) + l-(l + l) + 3x 1 + 1-1 



Indeed, in this case s± = 2 and Sj 2 = 1, hence there are 11 complex parameters in the soliton 
solutions (9 in the three \p) vectors, plus the two zeros k\ and k 2 ). The invariance matrix B can be 
found from the general formula ()91|) as 



B 



\ 



fell fei2 

fen 
fe 



\ 



J 



which has three free complex parameters. Thus A/f ree =11 — 3 = 8 as calculated above. 

Three solutions, with k x = 1 + i, k 2 = — 1 + 0.5i, {p^) — [ — 1, i, 1 — i] T ', and three different 
sets of |Poi ) an d |Poi ) vectors, are displayed in Fig. 3. In the generic case (first row in Fig. 3), 
this solution describes the breakup of a higher-order pumping wave (113) into the other u\ and u 2 
components (both higher-order). In non-generic cases, it can describe processes such as W2(sech) + 
«3 (higher-order) — > u\ (higher-order) + u 2 (higher-order) (second row of Fig. 3), Wi(sech) + W2(sech) + 
«3(sech) — > Ui (higher-order) + u 2 (higher-order) (last row of Fig. 3), and many others. The inverse 
processes of Fig. 3 can also be described by choosing cij and bj values such that u\ < u 3 < u 2 instead 
of u 2 < U3 < u\ in Fig. 3. 



V. CONCLUSION AND DISCUSSION 

We have proposed a unified and systematic approach to study the higher-order soliton solutions of 
nonlinear PDEs integrable by the iV x iV-dimensional Riemann-Hilbert problem. We have derived the 
complete solution to the Riemann-Hilbert problem with an arbitrary number of higher-order zeros, 
and characterized the discrete spectral data. As a result, the most general forms of higher-order multi- 
soliton solutions have been obtained in nonlinear PDEs integrable through the iV x iV-dimensional 
Riemann-Hilbert problem. In other words, the most general reflection-less (soliton) potentials in the 
iV-dimensional Zakharov-Shabat operators have been derived. The eigenfunctions associated with 
these reflection-less potentials are readily available from our soliton matrices. We have applied these 
general results to the three-wave interaction system, and new higher-order soliton and two-soliton 
solutions have been presented. These solutions reveal new processes such as M1+M2+M3 <-> M1+M2+W3. 
They also reproduce previous solitons in 0, 3, 2?| as special cases. Our results can be applied to 



derive higher-order multi-solitons in the NLS equation and the Manakov equations as well, but this 
is not pursued in this article. 
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The results obtained in this paper are significant from both physical and mathematical points 
of view. Physically, our results completely characterized higher-order solitons and multi-solitons in 
important physical systems such as the three-wave interaction equation, the NLS equation and the 
Manakov equations. These higher-order solitons can describe new physical processes such as those 
displayed in Figs. 1-3. If these integrable equations are perturbed (which is inevitable in a real- world 
problem), our higher-order solitons then become the starting point for the development of a soliton- 
perturbation theory which could determine what happens to these higher-order solitons under external 
or internal perturbations [3, From the mathematical point of view, our results completely 

characterized the discrete spectral data of higher-order zeros in a general iV-dimensional Riemann- 
Hilbert problem. These results will be useful for many purposes such as proving the completeness 
of eigenfunctions in a iV-dimensional Zakharov-Shabat spectral problem with arbitrary localized 
potentials. The difficulty of such a proof is caused by higher-order zeros. With our results, this 
difficulty can be hopefully removed. 

From a broader perspective, our results are closely related to many other physical and mathemat- 
ical problems. For instance, the lump solutions in the Kadomtsev-Petviashvili I eq uation are given 
by the higher-order poles of the time-dependent Schrodinger equation. In l2l| . lump solutions 



corresponding to certain special higher-order poles were derived, but the most general lump solutions 
still remain an open question. Note that the time- dependent Schrodinger equation is an infinite- 
dimensional system compared to our present iV-dimensional Riemann-Hilbert system. But the ideas 
used in this paper might be generalizable to the time- dependent Schrodinger equation as well. This 
remains to be seen. 
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APPENDIX A: GENERAL RIEMANN-HILBERT PROBLEM WITH ABNORMAL ZE- 
ROS 

Here we show that our soliton matrices of section IIHI can be generalized to the case of Riemann- 
Hilbert problem with abnormal zeros. However, due to the lack of important applications, we will 
show a simple example, which corresponds to a pair of zeros with different geometric multiplicities 
but the same algebraic multiplicity. Then we comment on the general case of several non-paired 
zeros. 

Let us use the simplest example to show the idea behind generalization of our results to the general 
Riemann-Hilbert problem with abnormal zeros. Consider one pair of zeros (ki,ki) which have the 
same algebraic multiplicity 2 but different geometric multiplicities which are 1 and 2 respectively. 
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The corresponding soliton matrices are as follows: 



k — ki 

r-\k) = (i+ (h^ihm) ( j + ^izSMM) , ( A2) 

with the conditions that (vj\vj) = 1, (t^l^i) = and (t>i|i>2) 7^ 0. To verify that the above matrices 
are indeed inverse to each other it is enough to rewrite the matrix T(k) in the form 

m = ( I+ (h-h)\v 2 )(v 2 \\ / 1 + (h-h^m^ 

V k — k\ J \ k — k\ J 

and take into account that Pj = \vj)(vj\ is a projector. Equations (jA2|) and ()A3|) are in fact the 
product representations of the form (j77|) . Now let us show that there are exactly two solutions to 
(p|r _1 (A;i) = 0. Indeed, the corresponding null vectors are as follows 

(Px\ = N, (p 2 | = N- (A4) 

This is due to the fact that r _1 (A; 1 ) = (I — Pi) (I — P2). But on the other hand there is just one 
solution to T(ki)\p) = 0: \pi) = \vi). Suppose that there is another solution \p 2 ) to T(ki)\p) = 
linearly independent from \pi). We have then using formula ()A1|) for T(ki): 

\P2) = \vi)(vi\p 2 ) + \v 2 )(v 2 \p2)- (A5) 

Thus IP2) = o\ v l) + ^l^)- Using this in formula ()A5|) we get, due to (u 2 |ui) = and (^i|t>2) 7^ 0, 

a\vi){vi\v 2 ) = 0, 

which is a contradiction, since a 7^ 0. 

The soliton matrices given by formulae (jAl|) - (jA2|) have the following form in the standard notations 
of Lemma ^ of section IIII1 _ _ _ _ 

m = I + \MM±^lM, ( A6) 

k — ki 

T-\k) = I + MM + MM + (AT) 
k — ki (k — kiy 

where 

|g"i) = (Aa - ^1)^2), |g 2 ) = (fci - fci)K), = (&i -^i) 2 (^iK)(^2|, (92! = (fci 



(fci - fci)(wi|tf 2 ) 

Notice that T(k) has two blocks of size 1, while T~ 1 (k) has one block of size 2. In general, for one pair 
of zeros with different geometric multiplicities, the soliton matrices have the structure of Lemma Q 
but with different numbers of blocks in T(k) and T~ 1 (k), while the total number of the \p)- and 
(p|-vectors appearing in these matrices is the same and equals to the order of the pair of zeros. One 
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can proceed to derive the representations similar to those in Lemma 0] for this case. Evidently, due 
to the way of the derivation, the formulae will be similar with the only difference in the number of 
blocks and block sizes in T(k) and r _1 (fc). 

In the more general case of the Riemann-Hilbert problem with abnormal zeros, the zeros can be 
non-paired (for instance, zero of order 2 in C + and two simple zeros in C_). Formally, this case can be 
obtained by "splitting" some of the zeros inside pairs in several distinct zeros in the soliton matrices 
r(fc) and T~ 1 (k) discussed above, since this limit is obviously regular (the geometric multiplicity of 
the zero to be split should be at least equal to the number of the generated in this way new zeros, 
thus providing for the needed number of blocks; formula (|A6|) . for instance, allows splitting of the 
zero k = k\ of T~ 1 {k) into two simple zeros). Thus, the most general case can be handled starting 
from the case of just one pair of zeros, i.e., the case discussed above. The explicit expressions for the 
soliton matrices T(k) and r _1 (A;) will involve similar relations between the numbers of zeros, their 
geometric multiplicities and the numbers and sizes of the //-blocks of vectors as those in lemma ^ 
though, obviously, with different particular numbers for each of the two matrices. 
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t=-15 



t=0 



t=15 




FIG. 1: Soliton solutions in the three- wave system (|14[) corresponding to a single pair of zeros with rank 
sequence {1, 2} at time t = —15,0 and 15. Here, k\ = 1 + i, \p^) = [~ 1> * 5 1 — i] T , \Pqi) = [1)0.5,— 1] T . 
First row: \p$) = [1, 1 + i, 0.5] T ; second row: \p$) = [0, 1 + i, 0.5] T ; third row: \p$) = [1,0,0.5] T . 
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t=-15 



t=0 



t=15 




FIG. 2: Soliton solutions in the three-wave system (|14[) corresponding to two pairs of simple zeros at time 
t = -15,0 and 15. Here, ki = 1 + i, k 2 = -1 + 0.51 First row: \p$) = [1, 1 + i, 0.5] T , |p£ } ) = [1,0.5, -1] T ; 
second row: \p$) = [0, 1 + », 0.5] T , = [1,0.5, -1] T ; third row: (p^) = [0, 1 + i, 0.5] T , IpjJ) = [1,0.5,0] T . 
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FIG. 3: Soliton solutions in the three-wave system (|14|) corresponding to two pairs of zeros — one elementary 
with algebraic multiplicity 2, and the other one simple. Here, k\ = 1 + i (elementary zero), &2 = — 1 + 0.5z 
(simple zero) , and \p^) = [ — 1, i, 1 — i] T ■ First row: \pqi) = [1, 1 + i, 0.5] T , \p^i) = [1,0.5,— 1] T ; second row: 
\p$) = [0,l + i,0.5] r , |p<?) = [1,0.5,-1] T ; third row: \p$) = [0, 1 + i, 0.5] T , |jsg? > = [1,0.5,0] T . 
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